Clustering, Chaos and Crisis in a Bailout Embedding Map 
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Abstract 

' We study the dynamics of inertial particles in two dimensional incompressible flows. The particle 

o ; 

■ dynamics is modelled by four dimensional dissipative bailout embedding maps of the base flow 



which is represented by 2 — d area preserving maps. The phase diagram of the embedded map 
is rich and interesting both in the aerosol regime, where the density of the particle is larger than 
that of the base flow, as well as the bubble regime, where the particle density is less than that 



Q . of the base flow. The embedding map shows three types of dynamic behaviour, periodic orbits, 



chaotic structures and mixed regions. Thus, the embedding map can target periodic orbits as well 
as chaotic structures in both the aerosol and bubble regimes at certain values of the dissipation 
parameter. The bifurcation diagram of the 4 — d map is useful for the identification of regimes 
^ , where such structures can be found. An attractor merging and widening crisis is seen for a special 

I region for the aerosols. At the crisis, two period-10 attractors merge and widen simultaneously 

■ into a single chaotic attractor. Crisis induced intermittency is seen at some points in the phase 

. diagram. The characteristic times before bursts at the crisis show power law behaviour as functions 

f-*) I of the dissipation parameter. Although the bifurcation diagram for the bubbles looks similar to 



that of aerosols, no such crisis regime is seen for the bubbles. Our results can have implications 
for the dynamics of impurities in diverse application contexts. 
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I. INTRODUCTION 



Inertial particle dynamics, i.e. the dynamics of small spherical particles immersed in fluid 
flows, has attracted considerable attention in recent years These studies 

have two fold importance. From the point of view of fundamental physics, the motion is 
governed by dynamical equations which exhibit rich and complex behaviour. From the point 
of view of applications, this dynamics constitutes the simplest model for impurities whose 
transport in flows is of practical interest, in contexts as varied as the behaviour of aerosols 
or pollutants in air, with consequences for weather and climate, and that of nutrients such 
as plankton in the sea with consequences for ocean life. 

The Lagrangian dynamics of such small spherical tracers in two dimensional incompress- 
ible fluid flows is described by the Maxey-Riley equations [3]. These are further simplified 
under various approximations to give a set of minimal equations called the embedding equa- 
tions where the fluid flow dynamics is embedded in a larger set of equations which include 
the the differences between the particle and fluid velocities l|, yl]. Although the Lagrangian 
dynamics of the underlying fluid flow is incompressible, the particle motion is compressible 
Sj, and has regions of contraction and expansion. The density grows in the former giving 
rise to clusters and falls in the latter giving rise to voids. The properties of the base flow 
have important consequences for the transport and mixing of particles. Map analogs of 
the embedding equations have also been constructed for cases where the fluid dynamics is 
modelled by area-preserving maps which essentially retain the qualitative features of the 
flow j9,[l^. The embedded dynamics in both cases is dissipative in nature. 

Further complexity is added to the dynamics by the density difference between the parti- 
cles and the fluid. In the case of two-dimensional chaotic flows, it has been observed earlier 
that particles with density higher than the base flow, the aerosols, tend to migrate away 
from the KAM islands, while particles lighter than the fluid, the bubbles, display the op- 
posite tendency pjj. Neutrally buoyant particles also showed a similar result, wherein the 
particles settled into the KAM islands. Our study indicates that in addition to the density 
difference, the dissipation parameter of the system also has an important role to play in the 
dynamic behaviour of the aerosols and the bubbles. 

In this paper, we study the dynamics of passive, flnite size, inertial particles in underlying 
two dimensional incompressible flows. The dynamics is modelled by the four dimensional 
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dissipative embeddings of two dimensional area-preserving maps. We obtain the phase 
diagram of the system in the a — 7 space where a is the mass ratio parameter, and 7 is 
the dissipation parameter. The phase diagram shows rich structure in the a < 1 aerosol 
regime, as well as the a > 1 bubble regime. Unlike the earlier results mentioned above, that 
the bubbles tended to form structures in the KAM islands[ll] and the aerosols were pushed 
away we found that structures form for both bubbles and aerosols in certain parameter 
regimes due to the role of the dissipation parameter. These regimes can be identified from 
the bifurcation diagram of the four dimensional embedding map. Both the aerosol and 
bubble regimes in the phase diagram show regions where periodic orbits as well as chaotic 
structures can be seen in the phase space plots. In addition to these, fully or partially mixed 
regimes can also be seen in the both aerosol and bubble cases. Thus the dynamic behaviour 
of the inertial particles can be of three major types. 

The bifurcation diagrams of the embedding map show further interesting structure. The 
bifurcation diagrams of the aerosol case show a regime where an attractor-merging and 
widening crisis can be seen. At the crisis, two period-10 attractors merge and widen simul- 
taneously into a single chaotic attractor. The parameter values at which the crisis occurs, 
are identified using the Lyapunov exponent and the bifurcation diagram. The signature of 
the attractor widening crisis can be seen in the intermittency of the time series of the map 
variables. The characteristic time before bursts r vs 7c — 7 follows a power-law, r ~ (7c — 7)'' 
[1^ . where 7c is the critical value of the dissipation parameter 7 at which the crisis occurs. 
Thus, the clustering of advected aerosol particles is dependent on the dissipation scale. How- 
ever, unlike the case of the aerosols, the bubble region did not show any of the characteristic 
signatures of crisis, nor was any intermittency seen in the bubble region. 

Our results can have implications for practical problems such as the dispersion of pollu- 
tants by atmospheric flows, and catalytic chemical reactions. 

II. THE EMBEDDING MAP 

The equation of motion describing the dynamics of passive neutrally buoyant particles of 
finite size in flows is obtained in This equation is called the bailout embedding equation, 
and is given by. 
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|[i-u(x)] = -(A/ + Vu).[i-u(x)]. (1) 

The velocity of the particle x and the local fluid velocity u(x) match, when the velocity 
gradient Vu exceeds the dissipation A, yielding a positive (A/ + Vu). When (A/ + Vu) 
becomes negative, the particle will get detached from the local fluid trajectory and would 
not follow the fluid. Therefore, the difference between the particle and the sorrounding fluid 
velocity can exponentially grow or get damped depending on the factor (A/ + Vu). 

The map analog of the flow bailout in Eqn. [T] is obtained in Ref ^ as. 



where the base map is an area preserving map given by Xn+\ = T{xn)- This is the bailout 
embedding map. We note that, the function (A/ + Vu) in the flow takes the form of 
K{xn) = e~^dxT in the map case. 

For the inertial particle case, the density of the particle and the fluid do not match. The 
dynamics of inertial particles in a fluid flow is described by the equation 
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Here, the velocity of the fluid parcel is u(x, y, t), and that of the particle advected by the 
fluid flow is V = d'iL/dt. The quantity a is the mass ratio parameter a = ^Pf/{pf + 2pp), 
where p/ is the fluid density and Pp is the particle density. Thus the a < 1 situation 
corresponds to aerosols and the a > 1 situation corresponds to bubbles. The dissipation 
parameter 7 is related to a by the relation 7 = 2a/3St, where St is the Stokes number. The 
dissipation parameter 7 gives a measure of the expansion or contraction in the phase space 
of the particle' s dy namics. 

Map analog 17 1 of the flow described by Eq. [3] has been choosen in Ref. |l| to be , 



x„+2 - M(x„,+i) = e ^(ax„+i - M(x„)) (4) 

where the base map is represented by an area preserving map Xn+i = M(x„+i). This can 
be rewritten, 

X„+i = M(X„) + Sn 

6n+i = e"'^[ax„+i - M(x^)]. (5) 
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FIG. 1: (a) The phase space plot of the standard map for K = 2.0. (b) The plot for the cor- 
responding embedding map for aerosols (a = 0.8, 7 = 0.8) and (c) bubbles (a = 1.2,7 = 0-4) ■ 
Random intial conditions distributed uniformly in the phase space were chosen. 



The new variable 5 defines the detachment of the particle from the local fluid parcel. 
This is the bailout embedding map. When the detachment measured by the 5 is nonzero, 
the particle is said to have bailed out of the fluid trajectory. In the limit 7 ^ cxd and 
q; = 1, 5 — >^ and the fluid dynamics is recovered. So the fluid dynamics is embedded in the 
particle's equation and is recovered under appropriate limit. This map is dissipative with 
the phase space contraction rate to be e"^". 
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A. The Embedded Standard Map 



We choose the standard map [IS*], to be our base map as it is a prototypical area- 
preserving system is widely studied in a variety of problems of both theoretical and ex- 
perimental interest. 

The standard map also serves as the usual test bed for the study of various transport 



phenomena and their quantifiers IJ] . The study of impurity dynamics using the standard 
map as the base flow, can thus give us a handle on the way in which these phenomena affect 
the impurity dynamics. Earlier studies of impurity dynamics using embeddings also use the 
standard map as the base map 

The parameter K controls the chaoticity of the map. The chaoticity of the map is 
dependent on the initial condition and the map nonlinearity parameter K. The map equation 
is given by, 

x„+i = x„ + y„+i {Mod 1) 

K 

Yn+i = Yn + TT- sin(27rx„) {Mod 1) (6) 

ZTT 

This map is taken to be the base map M(x^) in Eq. (3). The phase space of this base 
standard map for K = 2.0 is plotted in Fig. [T] in the region x G [—1, l],y G [—1, 1]. The 
particle dynamics is governed by the 4 dimensional map , 

K 

Xn+i = Xn + — sin(27ry„) + 6^ 
Zn 

K 

yn+i = Xn + yn + — sin(27ry„) + 6^ 
Zn 

^l+i = e-^[ayn+i - {yn+i - S^)] (7) 

It is clear that this 4-dimensional map is invertible and dissipative. The embedding map has 
3 parameters /T, 0,7. We plot the phase space portrait of the system evolved with random 
initial conditions at i^' = 2.0, 7 = 0.5, and two values of a, one in the aerosol regime 
< a < 1 (Fig. [Hb)), and one in the bubble regime 1 < a < 3 (Fig. [U^c)). The phase space 
plot of the standard map at the same value of K in Fig. [U^a) shows islands and chaotic 
regions separated by invariant tori which constitute barriers to transport. We observe from 
Figs. [U^b), and [D^c) that both the aerosols and bubbles have broken the barrier posed by 
the invariant curve and have targetted the islands forming structures. Thus, it is clear that 
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FIG. 2: Typical bifurcation diagrams for the embedded standard map for K = 2.0. a) aerosols 
(a = 0.8) b) bubbles (a = 1.2) for dissipation parameter ranging from 7 = 0.0 to 1.0. 

the dissipation parameter 7 plays a vital role in the formation of structures. Therefore in 
order to identify the regimes in which clustering or mixing can take place, it is necessary 
to study the full phase diagram. The bifurcation diagram of the embedding map can also 
provide insights into the regimes that can be expected in the phase diagram. 

III. THE BIFURCATION DIAGRAM 

We plot the bifurcation diagram of the embedded standard map (Fig. [2]) for the parameter 
values K = 2.0 and 7 values ranging from to 1, and for two a values, a = 0.8 and a = 1.2 
which correspond to the aerosol and bubble cases respectively. The bifurcation diagrams 
show many regions of periodic and chaotic behaviour for both the bubble and aerosol cases. 
This behaviour is reflected in the phase space plots of the map. The phase space plots for 
the aerosol case with 7 = 0.8, 7 = 0.3, and the bubble case with 7 = 0.1 can be seen in Fig. 
[3l Regimes with periodic orbits, structures and mixed regimes can clearly be seen here. 

Thus the phase space plots show different types of dynamical behaviour depending on the 
7 and a values. The phase diagram of the system in a — 7 space is essential for identifying 
the parameter regimes where different kinds of dynamic behaviours can be seen. 
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FIG. 3: Phase space plots for fixed K = 2.0 showing, (a) a chaotic structure in the phase space 
for the aerosol case a = 0.9, 7 = 0.75; (b) a chaotic structure in phase space for the bubble case 
a = 1.05, 7 = 0.75; (c) the mixing region in the bubble case a = 1.7, 7 = 0.3. 

IV. THE PHASE DIAGRAM 

The phase diagram of the system can be constructed using the largest Lyapunov exponent 
as the characteriser and is shown in Fig. HI For the calculation of the Lyapunov exponent, 
the first 5000 iterates are taken as transients and the next 1000 iterates are stored, for 100 
random initial conditions, uniformly spread in the phase space. 

The a values range from the aerosol regime < a < 1, to the bubble regime 1 < a < 3. 
Tongue-like structures are seen in the parameter space. Three types of regimes are seen 
in the phase diagram (i) regimes with periodic structures (P), where the largest Lyapunov 
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FIG. 4: (Color Online) The phase diagram. The regimes with periodic structures are marked with 
P (blue). The regimes with chaotic structures are marked with C (green for aerosols and yellow for 
bubbles). The regimes with structures in a mixing background are marked with a '*'(light blue). 
The regimes with full mixing are marked with M (red). 

exponent is negative, (ii) regimes with chaotic structures (C), (iii) mixing regimes (M), 
and (iv) chaotic structures in a mixing background. The Lyapunov exponent takes positive 
values in regimes (ii), (iii) and (iv). In the aerosol regime and the bubble regime, periodic 
orbits like those of Fig. [U^b) and (c) are seen inside the tongues marked by P (blue - 
(online)), and chaotic structures like those in Fig. [3t^b) are seen in the regions marked by 
C (green -(online), in the aerosol regime and yellow - (online), in the bubble regime) of the 
phase diagram. The fully mixing regimes marked with M (red - (online), which exist in both 
aerosol and bubble regimes, are identified by the following procedure. At a given parameter 
value, the phase space is covered by a 100 x 100 mesh, and the number of boxes accessed 
by the iterates is counted. The initial conditions and the transients are same as those that 
were used for calculating the Lyapunov exponent. The regions of parameter space where the 
iterates access more than 99% of the phase space grid turn out to be fully mixing regions 
with no remanent of any chaotic structure seen in the phase space plot, and have been 
marked with an M. 

The regions where the iterates access > 90% and < 99% of the boxes in the phase 
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FIG. 5: (a) A typical phase space plot for the regions with occupancy < 99% and > 90%. Chaotic 
structures are seen in a mixing background. These regions are identified in the phase diagram - 
'*' (light blue -online). (b) A typical phase space plot with occupancy < 90% and 50% 60%. They 
correspond to the unshaded (blank) regions in the phase diagram. 

space have been marked with an asterisk (*) (light blue- (online)). The phase space plots 
corresponding to this region have a chaotic structure in a well mixed background (Fig. 
Et^a)). The tiny patches of unshaded (blank) regions in the phase diagram are regions where 
^max > and the trajectories access less than 90% of the grid in the phase space. Typically, 
the trajectories cover 50 — 60% of the grid in the phase space in these regions. Fig. [5]^b) 
shows a typical phase space plot for this case. Here, the iterates wander densely along the 
separatrices in the phase space, but large voids are seen in the phase space. However many 
iterates stick in the neighbourhood of the tiny dense structures seen in the middle of these 
voids. 



V. THE AEROSOL CASE: A CRISIS IN THE EMBEDDED STANDARD MAP 



The bifurcation diagram of the embedded standard map shows regimes where the at- 
tractor undergoes a sudden discontinuous change of type as the system parameter is varied. 



Such sudden discontinuous changes signal the onset of a crisis [12j. The usual crises seen in 
dissipative dynamics are of the following types. In the attractor destruction type of crisis, 
a chaotic attractor is destroyed as the parameter passes through its critical crisis value. In 
the attractor widening or interior crisis, the size of the chaotic attractor increases suddenly. 
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In the attractor merging crisis, two or more chaotic attractors merge to form one chaotic 
attractor. It is to be noted that the inverse of the above processes (i.e, the sudden creation, 
shrinking, or sphtting of a chaotic attractor) occur as the parameter is varied in the other 
direction. In the case of our embedding map, the location and type of the crisis can be easily 
identified from the bifurcation diagram (see Fig. [2]^b)). 

The most prominent change in the size of the attractor is seen near 7 = 0.40. Here, as 
the parameter 7 decreases from 0.41 to 0.40 the attractor widens suddenly. The pre-crisis 
attractor has orbits of finite period, whereas after the crisis, the trajectories access the full 
range in y. The phase space plots at 7 = 0.41 and 7 = 0.40, corresponding to the pre-crisis 
and post-crisis situation respectively, are seen in Fig. [HI In the pre-crisis situation there 
are two period-10 attractors, and trajectories are confined to one or the other attractor, 
depending on the initial condition of the trajectory. After the crisis, the two attractors 
merge and widen into one. Thus we see the occurence of an attractor merging and widening 
crisis. 

The exact parameter value at which the crisis occurs can be identified using standard 
quantifiers like the Lyapunov exponents, and also the bifurcation diagram. 

The close-up of the bifurcation diagram in Fig. [2](b) near 7 = 0.40 is shown in Fig. 
ini^a). As the dissipation parameter 7 is reduced from 7 = 0.41 through 7 = 0.40, a sudden 
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FIG. 6: The identifiers of crisis: (a) the bifurcation diagram (b) the largest Lyapunov exponent 
plot, both obtained in the vicinity of the crisis, plotted as a function of the dissipation parameter 
7 show a sudden change at the crisis where the dissipation parameter value jc = 0.40452 . 
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FIG. 7: The time series plots for a single trajectory for (a)the pre-crisis attractor, 7 = 0.41; (b) 
the post-crisis attractor, 7 = 0.40. 

widening is seen. The largest Lyapunov exponent, can also be used to identify the parameter 
value at crisis. Fig. Mjo) shows the Lyapunov exponent plotted in the vicinity of the crisis. 
The exponent clearly shows a sudden jump at 7 = 0.40452, which signifies the occurence 
of the crisis [15]. From the above identifiers, it is seen that, the critical parameter value at 
crisis is 7c = 0.40452 (to our numerical accuracy). 

A. Crisis induced intermittency 

We see the phenomenon of intermittency in the time series of the phase space co-ordinates 
in the neighbourhood of the attractor widening crisis. Let 7c be the parameter value where 
the attractor widening crisis occurs. In the pre-crisis region, (7 > 7c), trajectories with 
arbitrary initial conditions in the phase space hop between the two attractors for small 
transient times. After this brief period, they settle into one of the attractors for arbitrarily 
long times (Fig. H^a)). 

In the post-crisis region, (7 < 7c), trajectories starting from arbitrary initial conditions 
hop intermittently between the regions corresponding to the two pre-crisis attractors for 
arbitrarily long times and never settle into either of the attractors. The time interval between 
the hopping is random as can be seen from the Figure. [7](b). In the attractor merging and 
widening crisis, an orbit initially confined to only one of the pre-crisis attractors (Figure. 
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FIG. 8: The phase space plots (a) the pre-crisis attractor with two period-10 orbits; 7 
(b)the post-crisis attractor; 7 = 0.40. 
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El^a)), say A or B, is able to access the widened single attractor composed of the two pre-crisis 
attractors A and B (Figure. EJ^b)) after 7 is reduced to 7c. 

When a trajectory leaves the pre-crisis attractor and accesses the post-crisis attractor, 
a burst is said to have occurred. The trajectory settles back into the pre-crisis attractor 
after a time interval and a new burst is initiated after some time. The characteristic time 
for the orbit to stay in the pre-crisis attractor before a burst occurs is defined as r. For 
each value of 7 the distribution of the laminar lengths r can thus be obtained. The average 
characteristic time r in the neighborhood of the 7c was seen to follow a power law 12|, 



1" ~ (7c - 7)^- 



(8) 



where the exponent was (3 = —0.35. Fig. [9] shows the log-log plot of r vs (7c — 7) which 
shows this scaling behaviour. 

We note that crisis induced intermittency is seen at several points on the left edge of the 
periodic tongue in the aerosol regime of the phase diagram. Power law scaling behaviour 
was seen at all these points. Unlike the aerosol case, no intermittency or crisis was observed 
anywhere in the bubble region. 

The reason for this difference in behaviour of the aerosols and the bubbles may lie in the 
structure of the phase space. The invariant tori which act as barriers to transport in the 
neighbourhodd of elliptic fixed points break up at higher values of K leaving behind leaky 
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barriers or cantori. In the case of the motion of inertial partcles, the elhptic islands and their 
neighbourhoods act as centrifuges, by pushing away the heavier aerosols and entrapping the 
lighter bubbles [l6|. At high values of the dissipation 7 , the dissipation counteracts the effect 
of the centrifugal force on the aerosols, so that the they are trapped in the neighbourhood 
of the elliptic islands, despite the centrifugal force that pushes them out. However, at lower 
values of dissipation e.g. for our map, at 7 < 7c, where 7c is the criticalparameter at which 
crisis occurs, the centrifugal force pushes out the aerosols through the leaky barries leading 
to attractor widening and intermittency. In the case of the bubbles, both the dissipation 
and the centrifugal force work to trap the bubbles, and no crisis or intermittency is seen. 

As mentioned above, the standard map which has been studied here is the prototypical 
area-preserving map. Hence, we expect that the key result, viz. that the dissipation of the 
fluid as well as the buoyancy of the particles affect the formation of coherent structures and 
the existence of mixing regions will apply to generic area preserving system. However, the 
details of the map, such as the existence of Cantori, and the size of the resonance zones will 
affect details like the existence of phenomena like crisis, crisis induced intermittency etc. 



VI. CONCLUSION 



To summarize, the present work discusses the advection of passive scalars of finite size in 
an incompressible two-dimensional flow for situations where the particle density can differ 
from that of the fluid. The motion of the advected particles is represented by an embedding 
map with the area preserving standard map as the base map. 

The new embedded standard map is invertible and dissipative and contains three pa- 
rameters, the nonlinearity K, the dissipation parameter 7 and the mass ratio parameter a 
which measures the extent to which the density of the particles differs from that of the fluid. 
The phase diagram of the system in the a — 7 parameter space shows very rich structure. 
Three types of dynamical behaviours are seen periodic orbits, chaotic structures, and mixed 
regimes which can be partially or fully mixed. Periodic structures are seen inside tongues in 
the parameter space in the aerosol as well as bubble regimes. Thus, the clustering or prefer- 
ential concentration of inertial particles can be seen in both the aerosol and bubble regimes. 
Chaotic structures can be seen outside the tongues in the aerosol regime. The bifurcation 
diagram indicates the existence of crises at parameter values on the edges of the tongues in 
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FIG. 9: The log-log plot of the characteristic time r versus 7c — 7 with the slope {3 = —0.35 it 
0.065087. 

this regime. The crisis is of the attractor merging and widening type in which two period 
ten orbits merge and widen into a chaotic attractor. Crisis induced intermittency is seen in 
the aerosol regime with charactersitic times r which scale as (7c — 7)^ where the exponent 
f3 = —0.35, and 7c is the critical value at which the crisis occurs. A well mixed regime is 
also seen for the aerosols. The bubble regime of the phase diagram shows periodic orbits, 
chaotic structures as well as partially and fully mixed regimes. No crisis or intermittency is 
seen here. 

The obvious generalisation of our study is to the three dimensional case volume-preserving 
case. In this case, due to the presence of additional degrees of freedom, a much richer phase 
diagram is expected. While earlier studies show that regions which target tubular structures 
in phase space do exist the presence of additional stretching directions would perhaps 
indicate much richer structure and larger mixing regions in phase space. Detailed studies in 
this direction are in progress. 

Our results have implications for the preferential concentration of inertial particles in 
flows, as well as for the targetting of periodic structures. Earlier results indicated that 
bubbles could breach elliptic islands and target structures, whereas aerosols could not. Our 
results indicate that the dissipation parameter 7 plays as crucial role in this as the den- 
sity differential parameter a, and the examination of the full phase diagram is necessary 
to draw conclusions about the parameter regimes where targetting and breaching can oc- 
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cur. We also see that mixing regimes can also exist, where random initial conditions can 
spread throughout the phase space. These results can have implications for the dynamics 
of impurities in diverse application contexts, e.g. that of the dispersion of pollutants in the 
atmosphere, flows with suspended microstructures, coagulation of material particles in flows 
and catalytic chemical reactions. We hope the present work will prove to be useful in some 
of these contexts. 
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